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Abstract 

A nonlinear wave alternative for the standard Black-Scholes option-pricing model is 
presented. The adaptive-wave model, representing controlled Brownian behavior of fi- 
nancial markets, is formally defined by adaptive nonlinear Schrodinger (NLS) equations, 
defining the option-pricing wave function in terms of the stock price and time. The model 
includes two parameters: volatility (playing the role of dispersion frequency coefficient), 
which can be either fixed or stochastic, and adaptive market potential that depends on 
the interest rate. The wave function represents quantum probability amplitude, whose 
absolute square is probability density function. Four types of analytical solutions of 
the NLS equation are provided in terms of Jacobi elliptic functions, all starting from 
de Broglie's plane-wave packet associated with the free quantum-mechanical particle. 
The best agreement with the Black-Scholes model shows the adaptive shock-wave NLS- 
solution, which can be efficiently combined with adaptive solitary-wave NLS-solution. 
Adjustable 'weights' of the adaptive market-heat potential are estimated using either 
unsupervised Hebbian learning, or supervised Levenberg-Marquardt algorithm. In the 
case of stochastic volatility, it is itself represented by the wave function, so we come to 
the so-called Manakov system of two coupled NLS equations (that admits closed-form 
solutions), with the common adaptive market potential, which defines a bidirectional 
spatio-temporal associative memory. 

Keywords: Black-Scholes option pricing, adaptive nonlinear Schrodinger equation, 
market heat potential, controlled stochastic volatility, adaptive Manakov system, 
controlled Brownian behavior 
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1 Introduction 



The celebrated Black-Scholes partial differential equation (PDE) describes the time-evolution 
of the market value of a stock option [TJ [2]. Formally, for a function u = u(t, s) defined on 
the domain < s < oo, < t <T and describing the market value of a stock option with the 
stock (asset) price s, the Black-Scholes PDE can be written (using the physicist notation: 
d z u = du/dz) as a diffusion-type equation: 

d t u = -^(as) 2 d ss ii - rs d s ii + ru, (1) 

where a > is the standard deviation, or volatility of s, r is the short-term prevailing 
continuously-compounded risk-free interest rate, and T > is the time to maturity of the 
stock option. In this formulation it is assumed that the underlying (typically the stock) 
follows a geometric Brownian motion with 'drift' \i and volatility a, given by the stochastic 
differential equation (SDE) [5] 

ds(t) = ns(t)dt + as{t)dW(t), (2) 

where W is the standard Wiener process. 

The economic ideas behind the Black-Scholes option pricing theory translated to the 
stochastic methods and concepts are as follows (see [6]). First, the option price depends on 
the stock price and this is a random variable evolving with time. Second, the efficient market 
hypothesis [8], i.e., the market incorporates instantaneously any information concerning 
future market evolution, implies that the random term in the stochastic equation must be 
delta-correlated. That is: speculative prices are driven by white noise. It is known that any 
white noise can be written as a combination of the derivative of the Wiener process |39] and 
white shot noise (see [9]). In this framework, the Black-Scholes option pricing method was 
first based on the geometric Brownian motion [H [2] , and it was lately extended to include 
white shot noise. 

The Black-Sholes PDE (pQ) is usually derived from SDEs describing the geometric Brow- 
nian motion (|2|), with the stock-price solution given by: 

In mathematical finance, derivation is usually performed using Ito lemma [37] (assuming that 
the underlying asset obeys the Ito SDE), while in physics it is performed using Stratonovich 
interpretation (assuming that the underlying asset obeys the Stratonovich SDE [38]) [9j [6]. 

The PDE ([1]) resembles the backward Fokker-Planck equation (also known as the Kol- 
mogorov forward equation, in which the probabilities diffuse outwards as time moves for- 
wards) describes the time evolution of the probability density function p = p(t, x) for the 
position x of a particle, and can be generalized to other observables as well [3]. Its first use 
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was statistical description of Brownian motion of a particle in a fluid. Applied to the option- 
pricing process p = p(t, s) with drift D\ = D\(t, s), diffusion D 2 = D2(t, s) and volatility a 2 , 
the forward Fokker-Planck equation reads: 

d tP = ^d ss [D 2 a 2 p] - 8 S (D lP ) . 

The corresponding backward Fokker-Planck equation (which is probabilistic diffusion in re- 
verse, i.e., starting at the final forecasts, the probabilities diffuse outwards as time moves 
backwards) reads: 

dtp = --<? 2 d ss (D 2 p) - d s (Dip) . 

The solution of the PDE (pQ) depends on boundary conditions, subject to a number of 
interpretations, some requiring minor transformations of the basic BS equation or its solution. 

The basic equation ([T]) can be applied to a number of one-dimensional models of inter- 
pretations of prices given to u, e.g., puts or calls, and to s, e.g., stocks or futures, dividends, 
etc. The most important examples are European call and put options (see Figured]), defined 
by: 

u C aii(s, t) = sjV(di) e~ TS - kN{& 2 ) e" rT , (3) 
n Put (s, t) = kAf(-d 2 ) e~ rT - sJV(-d a ) e~ TS , (4) 



AA(A) = ~ ( 1 +cri: 



a 



where erf(A) is the (real-valued) error function, k denotes the strike price and 5 represents 
the dividend yield. In addition, for each of the call and put options, there are five Greeks 
(see, e.g. [22]), or sensitivities of the option-price with respect to the following quantities: 

1. The stock price - Delta: AcaJi = d s ucaii and Ap ut = d s up u t, 

2. The interest rate - Rho: Pcall = ^r-^Caii and Pp ut = d r up u t] 

3. The volatility: Vega Ca n = dsucaii and Vega Put = d a u Put ; 

4. The elapsed time since entering into the option - Theta: 
6caii = cancan and P^t = d T up ut ; and 

5. The second partial derivative of the option-price with respect to the stock price - 
Gamma: r Ca ii = d ss u Ca ii and r Put = d ss up ut . 
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Figure 1: European call ([3]) and put options, as the solutions of the Black-Sholes PDE 
CD). Used parameters are: a = 0.3, r = 0.05, k = 100, 5 = 0.04 [22]. 

In practice, the volatility a is the least known parameter in (pQ), and its estimation is 
generally the most important part of pricing options. Usually, the volatility is given in a 
yearly basis, baselined to some standard, e.g., 252 trading days per year, or 360 or 365 
calendar days. However, and especially after the 1987 crash, the geometric Brownian motion 
model and the BS formula were unable to reproduce the option price data of real markets. 

The Black-Scholes model assumes that the underlying volatility is constant over the 
life of the derivative, and unaffected by the changes in the price level of the underlying. 
However, this model cannot explain long-observed features of the implied volatility surface 
such as volatility smile and skew, which indicate that implied volatility does tend to vary with 
respect to strike price and expiration. By assuming that the volatility of the underlying price 
is a stochastic process itself, rather than a constant, it becomes possible to model derivatives 
more accurately. 

As an alternative, models of financial dynamics based on two-dimensional diffusion pro- 
cesses, known as stochastic volatility (SV) models [ID] , are being widely accepted as a rea- 
sonable explanation for many empirical observations collected under the name of 'stylized 
facts' [llj. In such models the volatility, that is, the standard deviation of returns, originally 
thought to be a constant, is a random process coupled with the return in a SDE of the form 
similar to ([2]) , so that they both form a two-dimensional diffusion process governed by a pair 
of Langevin equations [T0l [HI [13] . 

Using the standard Kolmogorov probability approach, instead of the market value of an op- 
tion given by the Black-Scholes equation ([T]), we could consider the corresponding probability 
density function (PDF) given by the backward Fokker-Planck equation (see |9j). Alterna- 
tively, we can obtain the same PDF (for the market value of a stock option), using the 
quantum-probability formalism |44[ 145] . as a solution to a time-dependent linear Schrodinger 
equation for the evolution of the complex- valued wave ^—function for which the absolute 
square, \ip\ 2 , is the PDF (see [38]). 
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In this paper, I will go a step further and propose a novel general quantum-probability 
based option-pricing model, which is both nonlinear (see [Ml [XT] [TH] [20] ) and adaptive (see 
[21\ \19\ HI HO]). More precisely, I propose a quantum neural computation [32] approach to 
option price modelling, based on the nonlinear Schrodinger (NLS) equation with adaptive 
parameters. 

2 Adaptive nonlinear Schrodinger equation model 

This new adaptive wave-form approach to financial modelling is motivated by: 

1. Modern adaptive markets hypothesis of A. Lo |49[ I50|; 

2. My adaptive path integral approach to human cognition 29, 30], 13 lj ; 

3. Elliott wave (fractal) market theory [51 1 1521 [53]; and 

4. My recent monograph: 'Quantum Neural Computation' [32j . as well as papers on en- 
tropic crowd modelling based on the concept of controlled Brownian motion [33] 134"! 135] . 

To satisfy both efficient and behavioral markets, as well as their essential nonlinear com- 
plexity, I propose an adaptive, wave-form, nonlinear and stochastic option-pricing model with 
stock price s, volatility a and interest rate r. The model is formally defined as a complex- 
valued, focusing (1+1)— NLS equation, defining the option-price wave function ip = tjj(s,t), 
whose absolute square \ip(s,t)\ 2 represents the probability density function (PDF) for the op- 
tion price in terms of the stock price and time. In natural quantum units, this NLS equation 
reads: 



where dispersion frequency coefficient a is the volatility (which can be either constant or 
stochastic process itself), while Landau coefficient f3 = f3(r, w) represents the adaptive market 
potential, which is, in the simplest nonadaptive scenario, equal to the interest rate r, while in 
the adaptive case depends on the set of adjustable parameters {wi}. In this case, f3(r, w) can 
be related to the market temperature (which obeys Boltzman distribution [E]). The term 
V(ip) = — f3\i/j\ 2 represents the ip— dependent potential field. Physically, the NLS equation 
([5]) describes a nonlinear wave-packet defined by the complex- valued wave function ip(s,t) of 
real space and time parameters. In the present context, the space-like variable s denotes the 
stock (asset) price. 

1 Note that the domain of validity of the 'quantum probability' is not restricted to the microscopic world 
[46] . There are macroscopic features of classically behaving systems, which cannot be explained without 
recourse to the quantum dynamics (see [47] and references therein). 





(5) 
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2.1 Analytical NLS-solution 




Figure 2: The general Jacobi sine solution (|12[) of the NLS equation ([5]) with k = 1.2, m = 0.5, 
cr = /3 = 1, for t € (0,5) and s S (—7, 18). Thick line represents +Re[ip(s,t)], while dash line 
represents -Re[^(s, t)]. 

In case of (3 <C 1, V^) — > 0, so equation ([5|) can be approximated by a linear wave 
packet, defined by a continuous superposition of de Broglie's plane waves, associated with a 
free quantum particle. This linear wave packet is defined by the simple linear Schrodinger 
equation with zero potential energy (in natural units): 

id t il> = -^d ss ip. (6) 

Thus, we consider the ip— function describing a single de Broglie's plane wave, with the wave 
number (or, momentum) k and circular frequency u: 

i)(s,t) = 4){i)e i(ks - wt) , with £ = s - akt and <f>(£) G R. (7) 

Its substitution into the linear Schrodinger equation ([6]) gives the linear harmonic oscillator 
ODE, whose eigenvalues are natural frequencies of (|6|) and the solution is given by a Fourier 
sine or cosine series (see, e.g. [57| [58]). 

Similarly, substituting ([7]) into the NLS equation (0), we obtain a nonlinear oscillator 
ODE: 

+ [w - \ok 2 \ m + /30 3 (O = 0. (8) 
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Figure 3: The dark shock-wave solution (fT3j) of the NLS equation <J5J) with k = 1.2, a = (3 = 1, 
for £ G (0,5) and s 6 (—7, 18). Thick line represents +Ke[tp(s,t)], while dash line represents 
-RehHM)]. 

Following [55], I suppose that a solution </>(£) for ([8]) can be obtained as a linear expansion 

0(£) = d + aisn(0, (9) 

where sn(s) = sn(s,m) are Jacobi elliptic sine functions with elliptic modulus m € [0,1], 
such that sn(s, 0) = sin(s) and sn(s, 1) = tanh(s)H Using standard identities with associated 
elliptic cosine functions cn(£) and elliptic functions of the third kind dn(£), we have 



</>'(£) = a lC n(0dn(e), 

= -ol{sn(0[l - m 2 sn 2 (e)] + m 2 sn(£)[l - sn 2 (£)]}. 



Substituting (|9|) and (|10p into ([S]) ; after doing some algebra, we get 
a = 0, 



"j-±»'\/^p w = ^(1 + m 2 + /c 2 ), 



(10) 



(11) 



2 For example, the general pendulum equation: 

a"(t,<f>) + sin[a(t, <f>)] = 



has the elliptic solution: 



-i 



t, sin 
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Figure 4: The +tanh solution from Figure O with stochastic volatility at (random walk). 



-a 



which, substituted into the nonlinear oscillator 



±™>\/ — sn(£), 



-a 







tanh(£), 



gives 

for m £ [0, 1]; and 
for m = 1. 



Using the substitutions ((7|) and (fTT|) . we now obtain the exact periodic solution of ([5]) as 
^l(«,t) = ± m J-^ sn(s - akt) e i[ks-^t(i+ m 2 +k 2 )] ? 

V> 2 (s, t) = ■ 



-a 



tanh(s — akt) e 



i[ks-\at(2+k 2 )] 



for m G [0,1); (12) 
for m = 1, (13) 



where (fT2|) defines the general solution (see Figure [2|) , while (fT3|) defines the envelope shock- 
wav^ (or, 'dark soliton') solution (Figure [3|) of the NLS equation ([5]). The same shock-wave 
solution with stochastic volatility at (defined as a simple random walk) is given in Figure [H 



3 A shock wave is a type of fast-propagating nonlinear disturbance that carries energy and can propagate 
through a medium (or, field). It is characterized by an abrupt, nearly discontinuous change in the character- 
istics of the medium. The energy of a shock wave dissipates relatively quickly with distance and its entropy 
increases. On the other hand, a soliton is a self-reinforcing nonlinear solitary wave packet that maintains its 
shape while it travels at constant speed. It is caused by a cancelation of nonlinear and dispersive effects in 
the medium (or, field). 
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Figure 5: The general Jacobi cosine solution (|14p of the NLS equation ([5]) with k = 1.2, m = 
0.5, cr = /3 = 1, for t G (0,10) and s £ (—7,18). Thick line represents +Re[^(s,i)], while 
dash line represents -Re[ip(s,t)]. 

Alternatively, if we seek a solution as a linear expansion of Jacobi elliptic cosine 
functions, such that cn(s,0) = cos(s) and cn(s, 1) = sech(s)|j in a linear form: 

4>(C) = ao + aicn(£), 

then we get 



ip 3 (s,t) = ±m^^-^cn{s-akt)e^ ks -^ t{1 - 2m2+k2) \ for m £ [0,1); (14) 



r/> 4 (s,i) = ±^p^secHs-akt)e [[ks -^ t{k2 - 1)] , for m = 1, (15) 

where (|14p defines the general solution (Figure [5]) , while (|15p defines the envelope solitary- 
wave (or, 'bright soliton') solution (Figure [6]) of the NLS equation ([5]). The same soliton 
solution with stochastic volatility a t (a simple random walk) is given in Figure HI 



4 A closely related solution of an anharmonic oscillator ODE: 

(t>"(s) + < j,(s) + <t> 3 (s) = 

is given by 
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Figure 6: The bright solitary-wave solution (|15p of the NLS equation ([5]) with k = 1.2, 
a = (3 = 1, for t E (0, 10) and s £ (—7, 18). Thick line represents +Re[ip(s,t)], while dash 
line represents -Re[ip(s,t)]. 

In all four solution expressions (fl~2j) . ([13]) . ([T4"|) and (fToT) . the adaptive potential /3(w) is 
yet to be calculated using either unsupervised Hebbian learning, or supervised Levenberg- 
Marquardt algorithm (see, e.g. [411 142j). In this way, the NLS equation ([5]) resembles the 
'quantum stochastic- filtering neural network' model of [23] [2"4"1 125] . While the authors of 
the prior quantum neural network performed only numerical (finite-difference) simulations 
of their model, this paper provides theoretical foundation (of both single NLS network and 
coupled NLS network) with closed- form analytical solutions. Any kind of numerical analysis 
can be easily performed using above closed-form solutions ip^Sjt), (i = 1, ...,4). 

2.2 Fitting the Black— Scholes model using adaptive NLS— PDF 

The adaptive NLS-PDFs of the shock-wave type (|13p can be used to fit the Black-Scholes 
call and put options. Specifically, I have used the spatial part of (fl"3"l) . 




(16) 



where the adaptive market-heat potential /3(r, w) is chosen as: 




(17) 
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Figure 7: The +sech soliton from Figured! with stochastic volatility at (random walk). 



The following parameter estimates where obtained using 100 iterations of the Levenberg- 
Marquardt algorithm. In case of the call option fit (see Figure [8]), n = 5, 
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24.8952, 
24.8951, 
24.895, 
-37.3927, 
-37.2757, 
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2 



II':, 



II''. 



4 



-78112.3, 
-78112.3, 
-78112.3, u;jj 
= -3108.08, 
= -3968.35, 



wl 



a 



call 
NLS 



-0.119341ct BS 



Lcall 
ft NLS 



-48178.3, 
-48178.3, 
-48178.3, 
= -1520633, 
= -159782. 

'call 



0.0156422/cbs, 2nEs = 15.6423T BS . 



In case of the put option fit (see Figure ([9])), n = 3, 



wl = 0.000222367 
w 2 a = -0.428113 
wl = 4.70615 



irr 



= 82032.8, 
439.148, w. 



= 63876.9, 
205780.0, 



it') 



a 



put 
NLS 



27.1558, w 6 c = 139805.0 



-0.003444O-BS, k 



put 
NLS 



3.10354A:bs, T$£ s = -3103.54T BS . 



As can be seen from Figure ([9]) there is a kink near s = 100. This kink, which is a natural 
characteristic of the spatial shock- wave (|16p , can be smoothed out by taking the sum of the 
spatial parts of the shock-wave NLS-solution (fl~3|) and the soliton NLS-solution (fT5j) as: 



— [d\ tanh(s 

P 



kta) + c?2 sech(s — kta)] 



(18) 
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Figure 8: Fitting the Black-Scholes call option with /3(w)-adaptive PDF of the shock-wave 
NLS-solution (fTHjl . 

In this case, using 100 iterations of the Levenberg-Marquardt algorithm, the following pa- 
rameter estimates where obtained: 

w \ = -0.00190885, w\ = 6798.78, w\ = 5329.46, 
w\ = 18.1757, w\ = 23.5253, u% = 18354.9, 
w \ = -71.7315, wl = 4.15999, dj = 12807.2, 
d x = 0.345078, d 2 = -12.3948. 

^nls = -0.247932 BS , k»% = 0.260764/cbs, = 260.764T BS . 

The adaptive NLS-based Greeks can now be defined, using (3 = r and above modified 
(a, k, t) values, by the following partial derivatives of the spatial part of the shock-wave so- 



Delta = d s c/){s) = 2y^f J |f|sech 2 (s - fctcr) | tanh(s - kta)\abs' ( A /^tanh(s - kta)) , 



lution (H3]): 





| tanh(s — £;icr)| {crabs" tanh(s — kta)) 

f sinh(2s - 2kta)abs' (^™ tanh(s - kta))}], 
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Figure 9: Fitting the Black-Scholes put option with f3(w) -adaptive PDF of the shock-wave 
NLS ip 2 ( s it) solution (|13p . Notice the kink near s = 100. 



(-f) S/2 J\f\ tanh(s-fcto-)| tanh(s-Hff)|abs'(^f tanh(s-Ho-)) 



Rho = d r 6(s) 



Theta = $0(a) = (-p) 3 ^ J| ? |sech 2 (s - fctcr)| tanh(s - kta)\&bs' (y^f tanh(s - fctcr)) , 

where abs'(z) denotes the partial derivative of the absolute value upon the corresponding 
variable z. 



2.3 Coupled adaptive NLS— system for volatility + option-price evolution 
modelling 

For the purpose of including a controlled stochastic volatility into the adaptive-wave model 
©, the full bidirectional quantum neural computation model for option price forecasting 
can be represented as a self-organized system of two coupled self-focusing NLS equations: 
one defining the option-price wave function ip = i/j(s,t) and the other defining the volatility 
wave function a = a(s,t). The two NLS equations are coupled so that the volatility a is a 
parameter in the option-price NLS, while the option-price tp is a parameter in the volatility 
NLS. In addition, both processes evolve in a common self-organizing market heat potential, 
so they effectively represent an adaptively controlled Brownian behavior of a hypothetical 
financial market. 



5 Controlled stochastic volatility here represents volatility evolving in a stochastic manner but within the 
controlled boundaries. 
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Figure 10: Smoothing out the kink in the put option fit, by combining the shock-wave solution 
with the soliton solution, as defined by (fTHj) . 



Formally, I here propose an adaptive, symmetrically coupled, volatility + option-pricing 
model (with interest rate r and Hebbian learning rate c), which represents a bidirectional 
spatio-temporal associative memory. The model is defined by the following coupled-NLS+Hebb 
system: 

Volatility NLS : id t a = -^d ss a - (3 (\a\ 2 + 1 </>| 2 ) a, (19) 
Option price NLS : id t ip = --d ss tp - P (k| 2 + iVf) (20) 

N 

with: (3(r,w) = r^~^Wjgj, and 

i=i 

Adaptation ODE : Wi = —Wi + c\a\gi\ip\. (21) 

In this coupled model, the cr-NLS (|19j) governs the (s, t)— evolution of stochastic volatil- 
ity, which plays the role of a nonlinear coefficient in (|20p : the ^-NLS (|20p defines the 
(s, t)— evolution of option price, which plays the role of a nonlinear coefficient in (|19p . The 
purpose of this coupling is to generate a leverage effect, i.e. stock volatility is (negatively) 
correlated to stock returns^ (see, e.g. [H]). The w— ODE (|2~Tj) defines the (a, tp)— coupling 
based continuous Hebbian learning with the learning rate c. The adaptive market-heat po- 
tential f3(r,w), previously defined by (fT7]) . is now generalized into a scalar product of the 
'synaptic weight' vector Wi and the Gaussian kernel vector gi, yet to be defined. 

6 The hypothesis that financial leverage can explain the leverage effect was first discussed by F. Black [36] . 
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The bidirectional associative memory model (|19p - (|20p - (|2ip effectively performs quantum 
neural computation [32], by giving a spatio-temporal and quantum generalization of Kosko's 
BAM family of neural networks |26U27j . In addition, the shock- wave and solitary- wave nature 
of the coupled NLS equations may describe brain-like effects frequently occurring in financial 
markets: volatility /price propagation, reflection and collision of shock and solitary waves (see 
[28]). 

The coupled NLS-system (Tl9]l -(|20l). without an embedded w— learning (i.e., for constant 
(3 = r - the interest rate), actually defines the well-known Manakov system, which was 
proven by S. Manakov in 1973 [59] to be completely integrable, by the existence of infinite 
number of involutive integrals of motion. It admits 'bright' and 'dark' soliton solutions. 
Manakov system has been used to describe the interaction between wave packets in dispersive 
conservative media, and also the interaction between orthogonally polarized components in 
nonlinear optical fibres (see, e.g. |61[ I62j and references therein). 

The simplest solution of (fT9|) - ([20|) . the so-called Manakov bright 2-soliton, has the form 
resembling that of (US} and (Figure ED (see [63l M, EU ESI Ell EEl ES] ) , defined by: 

ip sol (s,t) = 26csech(26(s + 4at))e- 2i(2a2 * +as - 2b2 * ) , (22) 

where i/? sol (s,i) = f tj ) ' ° = ^ Cl,C2 ^ T * s a un ^ vec ^ or suc h that |ci| 2 + | C2 1 2 = 1. Real- 
valued parameters a and b are some simple functions of (a,f3,k), which can be determined 
by either Hebbian learning of Levenberg-Marquardt algorithm. Also, shock-wave solutions 
similar to (|13p are derived in Appendix. We can argue that in some short-time financial 
situations, the adaptation effect (|2~T1) can be neglected, so our option-pricing model (fT9|) - (i20|) 
can be reduced to the Manakov 2-soliton model ([22]) . as depicted and explained in Figure 

More complex exact soliton solutions have been derived for the Manakov system (I19p ~ 
([20]) with different procedures (see Appendix, as well as [7U [721 [73] ) . For example, in [71] . 
using bright one-soliton solutions (of the type of (fT5|) ) of the system (fT9l) - (i20|h many physical 
phenomena, such as unstable birefringence property, soliton trapping and daughter wave 
('shadow') formation, were studied. Similarly, searching for modulation instabilities and 
homoclinic orbits was performed in [7i \ \76\ [77j. In particular, local bifurcations of 'wave and 
daughter waves' from single-component waves have been studied in various forms of coupled 
NLS-systems, including the Manakov system (see [78 J and references therein). Let us assume 
that a small volatility a— component bifurcates from a pure option-price tp— pulse. Thus, at 
the bifurcation point, the volatility component is infinitesimally small, while the option-price 
component is governed by the equation 

d ss ip - i> + = o, 

whose homoclinic soliton solution is 

i/j(s) = \/2sechs. (23) 
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Figure 11: Hypothetical market scenario including sample PDFs for volatility |<r| 2 and \^\ 2 of 
the Manakov 2-soliton (I22p . On the left, we observe the (s, t)— evolution of stochastic volatil- 
ity: we have a collision of two volatility component-solitons, S\(s,t) and S2(s,t), which join 
together into the resulting soliton S2(s, t), annihilating the Si(s, t) component in the process. 
On the right, we observe the (s, t)— evolution of option price: we have a collision of two op- 
tion component-solitons, S\(s,t) and S2(s,t), which pass through each other without much 
change, except at the collision point. Due to symmetry of the Manakov system, volatility 
and option price can exchange their roles. 

A necessary condition for a local bifurcation of a homoclinic soliton solution with a small- 
amplitude volatility component from the option-price pulse (|23p is that there is a non-trivial 
localized solution to the linearized problem of the a— component. This takes the form of a 
linear Schrodinger equation 

d ss a - oj 2 a + 2sech 2 s a = 0, (24) 

which can be solved exactly (see |80j), and for local bifurcation we require a — ► as |s| — > 
±oo. 

As a final remark, numerical solution of the adaptive Manakov system (I19p -(j20 j) ~ (|2ip . 
with any possible extensions, is quite straightforward, using the powerful numerical method 
of lines (see Appendix in [32]). Another possibility is Berger-Oliger adaptive mesh refinement 
when recursively numerically solving partial differential equations with wave- like solutions, 
using characteristic (double- null) grids (see [79] and reference therein). 

2.4 Hebbian learning dynamics: analytical solution 

Regarding the Hebbian learning (f2Tj) embedded into the Manakov system (fl~9]) - (f20jh suppose 
e.g. that we have N = 10 synaptic weights (in a single neural layer), with the learning rate 
c = 0.7. The zero-mean Gaussians are defined as: 

_ jL 

gi = e ^, (i = l,...,N), 
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where {<7j} are (— 1, +1)— random standard deviations. Using random initial conditions, we 
get (by Mathematica of Maple ODE-solvers) the following analytical solutions of the Hebbian 
learning ODEs: 



m(t) = 

w 2 (t) = 
w 3 (t) = 
w±(t) = 
w 5 (t) = 
w 6 (t) = 
w 7 (t) = 
w 8 (t) = 

Wg(t) = 
W 10 (t) 



136485 erf(0.686579(106069 t - 1)) \ip\ 2 + 0.912318 \ip\ 2 - 0.00675663] 
0.932205 erf(0.553239(16336 t - 1)) |V>| 2 + 0.527646 |^| 2 - 0.249822] , 
0.471627 erfi(0.477447(2. 19341 t + 1)) |V| 2 - 0.274787 |V>| 2 + 0.582548] 
0.52899 erfi(0.6535(117079 t + 1)) |^| 2 - 0.453506 |^| 2 + 0.0773187] , 
0.362728 erfi(0.324902(4.73659 t + 1)) \ip\ 2 - 0.137812 \tp\ 2 + 0.798481] 
0.523292 erfi(0.6177(131043 t + 1)) |V| 2 - 0.416953 |^| 2 - 0.288671] , 
0.692907 erf(0.454319(2.42241 t - 1)) |^| 2 + 0.332217 |V>| 2 + 0.761879] 
0.432141 erf(0.315332(5.02843 t - 1)) |^| 2 + 0.148814 |^| 2 - 0.33264] , 



0.530916 erfi(0.673709(11016 t + 1)) \ip\ 2 - 0.473963 



0.17079 



e _t [0.443395 erf(0.322143(4.81806 t - 1)) |-0| 2 + 0.155768 |^| 2 + 0.49451] 



where erf(s) denotes the real-valued error function, while erfi(s) denotes the imaginary error 
function defined as: erf(is)/i. 




Figure 12: Time plot of the quick adaptive potential term (3(w) (as it appears in 
^j(s,f), i = l, ...,4) for the sample value of ip(s,t) = 0.5. 

In this way, we get the alternative expression for adaptive market-heat potential: (3(w) = 
r Tld=i w i9ii with interest rate r (see Figure [T2"j) . Insertion of 0(w), including the product 
\cr(s, t)\ \ip(s, t)\ calculated at time t into any Manakov solutions, gives the recursive QNN 
dynamics ip(s, t + 1) for volatility and option-price forecasting at time t+1. For example, an 
instant snapshot of the adaptive bright sech-soliton tp^s, t) is given in Figure [T3l 



3 Conclusion 

I have proposed a nonlinear adaptive-wave alternative to the standard Black-Scholes option 
pricing model. The new model, philosophically founded on adaptive markets hypothesis [39], 
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Figure 13: A snapshot of the adaptive ±sech-soliton ip 4 (s,t) with stochastic volatility at and 
trained potential j3{w) calculated at a sample fixed time to. We can see that due to quick 
learning dynamics, the whole solution is now decaying much faster than in Figure [6) 

50j and Elliott wave market theory |51}l52j. describes adaptively controlled Brownian market 
behavior, formally defined by adaptive NLS-equation. Four types of analytical solutions 
of the NLS equation are provided in terms of Jacobi elliptic functions, all starting from 
de Broglie's plane waves associated with the free quantum-mechanical particle. The best 
agreement with the Black-Scholes model shows the adaptive shock-wave NLS-solution, which 
can be efficiently combined with adaptive solitary-wave NLS-solution. Adjustable 'weights' 
of the adaptive market potential are estimated using either unsupervised Hebbian learning, 
or supervised Levenberg-Marquardt algorithm. For the case of stochastic volatility, it is 
itself represented by the wave function, so we come to the integrable Manakov system of 
two coupled NLS equations with the common adaptive potential, defining a bidirectional 
spatio-temporal associative memory machine. 

As depicted in most Figures in this paper, the presented adaptive-wave model, both 
the single NLS-equation ([5]) and the coupled NLS-system (fT9|) - ([20]) - ([2T|) . which represents a 
bidirectional associative memory, is a spatio-temporal dynamical system of great nonlinear 
complexity (see |40|). much more complex then the Black-Scholes model. This makes the 
new wave model harder to analyze, but at the same time, its immense variety is potentially 
much closer to the real financial market complexity, especially at the time of economic crisis 
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abundant in shock-waves. 

Finally, close in spirit to the adaptive-wave model is the method of adaptive wavelets in 
modern signal processing (see |81| and references therein, as well as [32J for an overview), 
which could be used for various market dimensionality reduction, signals separation and 
denoising as well as optimization of discriminatory market information. 



4 Appendix: Manakov system 

Manakov's own method was based on the Lax pair representation^ Alternatively, for nor- 
malized value of the market-heat potential, (3 = r = 1, Manakov system allows solutions of 
the form: 

a(s, t) = <p( a ) e^, V(*, t) = <j>(s) i<\ (26) 

where <p, cj> are real-valued functions and w a ,w^, are positive wave parameters for volatility 
and option-price. Substituting (126ft into the Manakov equations we get the ODE-system [62] 



<ff'{s) = wMs)-[^(s) + ^(s)) ( p(s), (27) 

cf>"(s) = wl<f>(s)-[cj> 2 (s) + <p 2 (s)}<t>(s). (28) 

For w a = = w, equations (|27p ~ (|28p have a one-parameter family of symmetric single- 
humped soliton solutions (see the left part of Figure [T^|) given by 

ip(s) = ±(j)(s) = w sech(u) s), (29) 

as well as periodic solutions: 

ip{s) = Acos(Bs) and <£(s) = A sin(Bs), (30) 

where A = \/w 2 + B 2 (with B an arbitrary parameter). For < w < 1 there is also another, 
in general asymmetric, one-parameter family of solutions for each fixed w [62] 



ip(s) = y/2{\ - w 2 )cosh (ws)/k, (31) 
4>(s) = —w\/2(l — w 2 ) sinh(s — sq)/k, where 
K = cosh(s) cosh(u; s) — w sinh(s) sinh(u; s), 



7 The Manakov system (|19p - (|2Up has the following Lax pair [7U] representation: 

d x (j} = M(j) and d t <f> = B<j>, or d x B — d t M = [M, B], with (25) 

/ -iA Vi i>2 \ 
M(A) = -ip 1 iA and 
\ -i>2 iA / 
/ 2A 2 - l^xl 2 - \ip 2 \ 2 2itp 1 X-d x ip 1 2iip 2 \~d x 4> 2 \ 
B(A) = -i -21^-0^ -2X 2 + \i} 1 \ 2 Vi^ 2 

V -2iip;\-d X Tp* -2A 2 + |^ 2 | 2 / 
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in which cp is symmetric and <f> antisymmetric. 




Figure 14: Initial envelopes for volatility \<?q\ and option-price |^ | within the Manakov 
solitons: (left) bright compound soliton (|29p . (middle) dark compound soliton (|32p . and 
(right) compound kink-shaped soliton (|36|) . These initial envelopes can be used for numerical 
studies of the Manakov system and its various generalizations (modified and adapted from 

On the other hand, for negative values of the potential 0, the Manakov equations accept 
dark soliton solutions of the form [75] 

a(s, t) = ^(s, t) = k [tanh(fes) - i] e^'-s* 2 *) j ( 32 ) 

which are localized dips on a finite-amplitude background wave (see the middle part of Figure 
[T4"|) . In this very interesting case, volatility and option-price fields are coupled together, form- 
ing a dark compound soliton. Note that their respective relative amplitudes are controlled 
by the corresponding nonlinearities and frequency. For (3 = — 1 the Manakov equations alow 
also solutions of the form: 

a(s, t) = ip(s) e- iw °\ ip(s, t) = <j){s) e - ™*'. (33) 

Introducing (|33p into the Manakov equations, we get the ODE-system: 

p"(s) = [<p 2 (s) + ^(s)]<p(s)-w 2 Ms), (34) 
0"( 5 ) = [^W + ^W^-w^s), (35) 

which, for w a = = w, allow for kink-shaped localized soliton solutions (see the right part 
of Figure [Til) given by [75] 

<p(s) = ±0(s) = (w/y/2) tanh(w s/y/2), (36) 

as well as periodic solutions (|30]l . Inserting ([TS|) back into (f33l) gives the double-kink solution 
for the Manakov system: 

a(s, t) = ±{w/V2) tanh(w a/y/2) e~ iwH , ip{s, t) = ±(w/V2) tanh(w s/\/2)e"™ 2 *. (37) 
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